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ABSTRACT 

We present a study of the turbulent velocity fields in the Intra Cluster Medium of a sample of 
21 galaxy clusters simulated by the SPH-code Gadget2, using a new numerical scheme where 
the artificial viscosity is suppressed outside shocks. The turbulent motions in the ICM of our 
simulated clusters are detected with a novel method devised to better disentangle laminar bulk 
motions from chaotic ones. We focus on the scaling law between the turbulent energy content 
of the gas particles and the total mass, and find that the energy in the form of turbulence scales 
approximatively with the thermal energy of clusters. We follow the evolution with time of the 
scaling laws and discuss the physical origin of the observed trends. The simulated data are in 
agreement with independent semi-analytical calculations, and the combination between the 
two methods allows to constrain the scaling law over more than two decades in cluster mass. 
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1 INTRODUCTION 

A good deal of evidence from observational and numerical works 
is suggesting that a non negligible budget of turbulent motions is 
stored in the Intra Cluster Medium (ICM) of galaxy clusters. The 
observational evidence from the gas pressure map of the Coma 
cluster (Schuecker et al. 2004), the lack of resonant scatterings in 
Perseus (Churazov et al. 2004), and the many indications of non- 
thermal emissions possibily due to acceleration of particles into 
a turbulent magnetized medium (e.g. Brunetti et al.2004 and ref- 
erences therein for a review), are pointing towards the relevance 
of chaotic motions within the ICM. Moreover, Eulerian numerical 
simulations of merging clusters (e.g., Bryan & Norman 1998; Roet- 
tiger. Stone & Burns 1999; Ricker & Sarazin 2001) or of accretions 
of less massive structures as dwarf galaxies (e.g. Mayer et al. 2005) 
have provided good representations of the way in which turbulence 
may be injected into the ICM. Also theoretical calculations expect a 
non-negligible amount of turbulent motions in galaxy clusters (e.g. 
Cassano & Brunetti 2005, C&B05), possibly with a relevant am- 
plification of seed magnetic fields through kinetic dynamo-effects 
(e.g. Subramanian, Shukurov & Haugen 2005; EnBlin & Vogt 2005; 
Schekochihin & Cowley 2006). Even the claim against turbulence 
in the ICM of some observed galaxy clusters as Perseus (e.g. Fabian 
et al. 2003) actually rule out only scenarios of very strong turbu- 
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lence, i.e. with an energy density in turbulence not far from the 
thermal support. 

Unluckily, due to the abrupt breaking of the main instrument 
onboard of the Astro-E2 mission, the direct detection of turbulent 
fields through the broadening of iron lines profile (e.g. Inogamov & 
Sunyaev 2003) has to be postponed to the future. From the numer- 
ical viewpoint, Smoothed-Particles-Hydrodynamics (SPH) repre- 
sents a unique tool to investigate the physics of both the collisional 
(gas) and non collisional (dark matter - DM) mass components 
of galaxy clusters; moreover, the high dynamic coverage of SPH 
permits to study a large interval of cluster sizes, and allows one 
to follow the evolution of the most interesting features with cos- 
mic time. In this Letter we present the first results from a quest 
for turbulent velocity fields in a sample of well-tested SPH sim- 
ulated galaxy clusters (Sec.|2j; so far, this is the first time that a 
fully cosmological and highly resolved set of simulations is em- 
ployed to search and characterize ICM turbulence. To this end we 
developed a novel approach to the problem of disentangling bulk 
motions of gas particles from chaotic ones (Sec.|3}; in Dolag et al. 
2005 (hereafter Do05) we claimed that such method is more effi- 
cient than only removing the cluster bulk velocity as usually done 
in literature (e.g. Bryan & Norman 1998), allowing one to focus 
only on the most chaotic part of the ICM flow. In this Letter we 
focus on the scaling relation between the thermal energy content 
of simulated galaxy clusters and their kinetic energy in the form of 
turbulent motions (Sec.|4} and finally compare our findings with 
the expectations from semi-analytical calculations (Sec.|5}- 
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Figure 1. 2-dimensional maps for the central ~ (l.SMpc)^ region of the most massive cluster of our sample. A subcluster ten times less massive is being 
accreted (entering from the left). The streamlines on the left depict the velocity defined with respect to the center of mass velocity; on the right depict the 
local velocity residuals for two different resolutions of the algorithm: I ~ 32kpc (center) and I ~ 16kpc(rig/i0. The overlayed colour maps show the kinetic 
energy of gas particles for each of the velocity fields. 



2 THE SAMPLE OF SIMULATED CLUSTERS. 

The detailed information about our cluster sample can be found 
in Do05. It consists of 9 resimulations with 21 galaxy clusters 
and groups simulated with the tree N-body-SPH code Gadget2 
(Springel 2005), performed several times with different physical 
processes. The cluster regions were extracted from a dark matter 
only simulation with box of 479 /i~^Mpc on a side and in the con- 
text of a ACDM model with fio = 0.3, h = 0.7, crs = 0.9 
and fib = 0.04 (Yoshida, Sheth & Diaferio 2001). Adopting the 
'Zoomed Initial Conditions' technique (Tormen, Bouchet & White 
1997) the regions were re-simulated in order to achieve higher 
mass and force resolution. Particle mass for the resimulations is 
todm = 1.13 X 10^ h-^MQ and mgas = 1.7 x 10** h'^MQ-, the 
gravitational softening adopted is fti ~ 5 ft"^kpc (Plummer equiv- 
alent between z = 5 and z = and kept fixed in comoving units at 
higher redshits; this also marks the minimum value for smoothing 
the gas particles within the SPH algorithm). 

Here we used all the 21 haloes having virial masses in the 
range Mvir = 5.3 • lO" - 2.3 • IO^^Mq ■ h'^ found in the set of 
simulations. They are resolved by 8 ■ 10* to 4 • 10® gas and DM par- 
ticles respectively. Even if the whole sample allows one to study the 
role of plasma conductivity, cooling, star formation and feedback 
processes, we restricted so far to a non-radiative SPH subset where 
an improved recipe for the numerical viscosity of gas-particles is 
used. This recipe, which follows an idea of Monaghan & Morris 
(1997) manages to reduce the amount of numerical, and therefore 
un-phyisical dissipation of chaotic motions, ensuring at the same 
time a good treatment of shocked features; its use is mandatory for 
any study of the chaotic part of the ICM dynamics, since otherwise 
the standard SPH scheme greatly suppresses any chaotic motion at 
the smallest scales, even in absence of shocks. For a more detailed 
discussion about this method we address the reader to Do05. 



3 DETECTION OF TURBULENT MOTIONS. 

In order to characterize turbulent velocity fields in a fluid medium, 
the crucial point is to extract a pattern of velocity fluctuations from 
a complex distribution of velocities. The most used approach in the 
literature (e.g. Norman & Bryan 1998; Sunyaev, Norman & Bryan 
2003) is simply to use individual velocities of the gas particles after 



subtracting of the mean velocity, computed within a fixed volume. 
Although this method has been widely employed in many previ- 
ous works, which indeed found a remarkable level of turbulence 
within simulated ICM, it can turn out to be highly misleading in 
the case of substructure crossing the cluster volume. As shown in 
the left panel of Fig. Q the motion of a subcluster usually tracks 
laminar velocity patterns which can differ from the mean virial ve- 
locity. This provides a spurious contribution to the estimated turbu- 
lent energy, whereas the only correct contribution to consider is the 
chaotic velocity field at the interface layer with the resident ICM 
of the primary cluster and along the tail of the subcluster In order 
to improve on this we conceived an algorithm which disentangles 
the chaotic part of the flow in a more effective way. We defined a 
mean local velocity field Vcdi in a cell by interpolating the velocity 
and density of each gas particle onto a regular mesh using a Trian- 
gular Shape Cloud (TSC) window function. The dependence of the 
mesh-spacing to the final results is discussed in section|4| We then 
evaluated the local velocity fluctuations of each gas particle by: 

5vi ~ Vi - Vcoll, (1) 

where Vi is the 3-dimensional particle velocity, Vcdi is the local 
mean velocity field of the cell where each particle falls. This field 
obviously varies depending on the choice of grid size: the subtrac- 
tion of a local velocity field is expected to progressively filter out 
the contribution from laminar bulk-flows produced by the gravita- 
tional infall (which is certainly relevant, for any redshift and all dis- 
tances from the cluster center, e.g. Tormen, Moscardini Sl Yosida 
2004). FigureQgives an example of the velocity fields obtained by 
subtracting the mean bulk velocity {left panel) and after subtract- 
ing the local mean velocity found by our new algorithm for two 
representative resolutions of the TSC algorithm {middle and right 
panel). The laminar flow at the largest scales is filtered out by the 
TSC-kernel, while the most chaotic and swirling flows are high- 
lighted, even though some of the biggest vortices might also have 
been filtered out. We stress here that this method is just a zero-order 
approximation: an accurate spectral analysis for turbulent dynam- 
ics within the simulated ICM will require further improvements, 
since the ICM is likely characterized by multiple energy injections 
(from mergers and accretions) driven simultaneously from different 
scales. 
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Figure 2. Individual patlis for fours clusters of the sample, in the 
log Energy — log Mtot plane. The upper two panels show the evolution 
of the most and of the less massive cluster within our catalogue, whereas 
the lower two panels show the evolution of two clusters with a nearly equal 
final mass (~ 5 X IO^^Mq/i"^), but very different "relaxation" state: left 
panel is for the "relaxed" (i.e. f < 0.5 at z=0) cluster g914 while the right 
one is for the "perturbed" one, g8 b ^ 0.5 at z=0). 



4 SCALING LAWS FOR TURBULENT KINETIC 
ENERGY 

The main goal of this section is to investigate the scaling laws be- 
tween the mass (gas plus dark matter particles) of clusters/groups, 
Aftot, and the thermal, kinetic and turbulent energy of the ICM. 

Due to computational limitations we so far restricted our anal- 
ysis to a cubic region, centered onto the center of the cluster, of 
equivalent volume Vbox ~ (i?vir)^- This ensures that we consider 
in any case a number of gas particles ranging from several thou- 
sands to nearly one million. After the velocity decomposition is 
performed (section|3}, we evaluate the turbulent energy content as: 



-Etur = -^mg 



(2) 



where the sum is done over the module of the velocity fluctuation, 
Svi, of the gas particles. This calculation was repeated at three dif- 
ferent resolutions of the TSC-kernel used to define the local mean 
velocity field: I =16, 32 and 64 kpc. As discussed in Do05 a rea- 
sonable subtraction of the laminar pattern is obtained with a mesh 
spacing around ~ 30kpc (corresponding to a TSC equivalent length 
three times larger)^. 

The total kinetic and thermal energies were evaluated as: 



-Eth = -m„ 



BOX 



flTUp 



(3) 



^ Note that this filtering width is always larger than the smoothing length of 
gas particles within the volume region we consider, for each cluster and each 
redshift of observation. An adaptive sampling scheme would be expecially 
mandatory to catch turbulence in the cluster outskirts, due to the decrease 
of particle density. 



Table 1. Values for the slopes of the kinetic and turbulent scaling laws at 
zero redshift, for the whole sample of data and the "relaxed" subsample, 
with Icr en'ors. 



I 


Dtur (all) 


Dtur (relax) 


Dk (all) 


Dk (relax) 


16 kpc 


1.43 ±0.06 


1.63 ±0.04 


1.38 ±0.04 


1.72 ±0.03 


32 kpc 


1.49 ±0.04 


1.72 ±0.01 


1.38 ±0.04 


1.72 ±0.03 


64 kpc 


1.49 ±0.03 


1.73 ±0.05 


1.38 ±0.04 


1.72 ±0.03 



where fcs is Boltzmann's constant, Ti the gas particle temperature, 
^ = 0.59 the mean molecolar weight in AMU, nip the proton mass 
and /e = 0.58 the fraction of free electrons per molecule, assuming 
a primordial mixture of xh = 0.76, and 
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where the module of velocity, Vi, has been reduced to the center of 
mass velocity frame (as in Norman & Bryan 1998). 

In figure|2|we report the time evolution of four representative 
clusters in our sample in the i5TUR,TH,K-Aftot plane. The most 
"relaxed" structures (as the cluster g914, bottom left panel) present 
a fairly smooth evolution, whereas "perturbed" structures (as g8b 
and glf, right panels) show a more complex evolution with episodic 
jumps in turbulent and kinetic energies, and have a high ratio, 5, 
between kinetic (and turbulent) and total (thermal plus kinetic) en- 
ergy. This reflects the significant difference in the ratio between the 
kinetic and the potential energy of these clusters (e.g. Tormen et 
al. 1997), which is higher for the perturbed ones. 

Since our cluster sample is extracted from re-simulations cen- 
tered on 9 massive and fairly isolated clusters, smaller systems gen- 
erally correspond to structures about to be accreted by larger ones. 
As such, small systems are often perturbed, and this introduces a 
bias in the dynamical properties of the cluster population. This bias 
can however be alleviated by restricting our analysis only to the 
most "'relaxed" objects in our sample, as we will see below. 
In general, we find the following power law scaling between cluster 
energy (thermal, kinetic or turbulent) and cluster mass: 



IQISMq fe- 



es) 



with j = TH, K, TU R, and where Aj and Dj are the zeroth point 
and the slope of the correlations, respectively. 
We find that the scaling of thermal energy with mass is always con- 
sistent with that expected in the virial case, Dth ~ 5/3, while 
the values of Dk and Dtur slightly depend on the number of 
"perturbed" small systems included in the analysis. With all sys- 
tem included, the slope of the scaling between turbulent energy and 
cluster mass is flatter than that between thermal energy and mass by 
~ 0.2. As we remove more and more small perturbed systems, the 
turbulent slope steepens toward the thermal value. If we define the 
ratio ^ between turbulent and thermal energy for each system and 
redshift, we find that the flattening of the turbulent scaling with re- 
spect to the thermal scaling is statistically significant only if objects 
with 5 ^ 0.5 (nine at z=0) are included. 

The slopes of the scalings thus obtained are stable and do not 
depend on the value of the mesh-spacing, /, adopted for the sub- 
traction of the laminar motions (Fig.|3}, as the maximum difference 
does not exceed ADtur ~ 0.1 at any redshift. 

In addition, the dependence of the turbulent energy on the 
adopted mesh-spacing confirms the behaviour already presented 
in Figure 4 of Do05, with Etvu. <xl^^^. 
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Figure 3. Scaling laws at redshift z = for the 1 2 most relaxed clusters 
(5 < 0.5); the values of the slopes for the different relations are reported in 
the panel. For the sake of displaying, only the datapoints of the I = 32kpc 
grid turbulence are drawn. 
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Figure 4. Redshift evolution of slope(/e/i') and zero point (right) of the scal- 
ing law Eq.js), for the sample of objects with ^ < 0.5. 



Finally, Figure |4| shows the redshift evolution of the slopes, 
Dj, and of the zero points, Aj, of the five correlations (Eq. |5}- 
It is clear that the slopes are relatively constant with redshift; this 
does not change significantly, unless very "perturbed" groups with 
^ ^ 0.5 (at each z) are considered in the analysis. In this last case 
a sistematic flattening (ADk,tur ~ 0.2) of the scaling of the 
kinetic and turbulent energies with cluster mass at low redshift is 
found; this is caused by the interactions between objects, which 
makes the smaller systems more and more perturbed as time pro- 
ceeds. 



5 COMPARISON WITH SEMI-ANALYTICAL RESULTS 

In the previous section we reported on the scaling between the tur- 
bulent energy and the thermal (and kinetic) energy as measured in 
simulated clusters, without motivating their physical origin. Clus- 
ter mergers are likely to be responsible of most of the injection 



of turbulent velocity fields in the ICM. Simulations have indeed 
shown that the passage of a massive sub-clump through a galaxy 
cluster can trigger turbulent velocity fields in the ICM (Roettiger et 
al.l999; Ricker & Sarazin 2001; Do05). 

A simple method to follow the injection of merger-turbulence 
during cluster life is given by semi-analytical calculations. C&B05 
used merger trees to follow the merger history of a synthetic popu- 
lation of galaxy clusters (using the Press & Schechter 1974 model) 
and calculated the energy of the turbulence injected in the ICM 
during the mergers experienced by each cluster. In these calcula- 
tions turbulence is injected in the cluster volume swept by the sub- 
clusters, which is bound by the effect of the ram pressure stripping, 
and the turbulent energy is calculated as a fraction of the PdV work 
done by the sub-clusters infalling onto the main cluster. 

Although simplified, this semi-analytical approach allows a 
simple and physical understanding of the scaling laws reported 
in the previous Section. Indeed, since the infalling sub-clusters 
are driven by the gravitational potential, the velocity of the infall 
should be ~ 1.5 — 2 times the sound speed of the main cluster; 
consequently, the energy density of the turbulence injected during 
the cluster-crossing should be proportional to the thermal energy 
density of the main cluster. In addition, the fraction of the vol- 
ume of the main cluster in which turbulence is injected (the vol- 
ume swept by the infalling subclusters) depends only on the mass 
ratio of the two merging clusters, provided that the distribution of 
the accreted mass-fraction does not strongly depend on the clus- 
ter mass (Lacey & Cole 1993). The combination of this two items 
yields a self-similarity in the injection of turbulence in the ICM; 
the energy of such turbulence should scale with the cluster thermal 
energy and the turbulent energy should scale with virial mass with 
a slope Ds^m ~ 1.67 (C&B05). 

In Fig|5|we report the integral of the turbulent energy (in- 
jected in the ICM up to the present time) versus the cluster mass, 
as estimated under the C&B05 approach with 360 merging trees 
of massive galaxy clusters, together with the measures done on 
our hydrodynamically simulated clusters; the two scalings are con- 
sistent within la errors. The two approaches are complementary, 
since semi-analytical calculations can follow the properties of > 
IO^^Mq clusters which are rare in numerical simulations due to 
the limited simulated cosmic volume, and so provide a physical ex- 
trapolation of the trend derived by our numerical simulations. This 
strengthens our claim that the turbulent velocity fields detected in 
simulated clusters are actually real turbulent fields supplied by the 
mass accretion process acting in galaxy clusters. 

Both estimates of the turbulent energy in galaxy clusters show 
an overall content of turbulence which ranges from 25% to 35% of 
the thermal one, in the {RvirY^ region. This should be considered 
as an upper limit of the turbulent energy content at a given time, 
because simulations do not contain appropriate recipes for the dis- 
sipation of the turbulent eddies at the smallest scales (for the sake 
of comparison, the semi-analytical calculations in Fig|5|were con- 
ceived to focus on the injection of turbulence in the whole cluster 
life). 

Fig|5| highlights the different behaviour of "perturbed" (i.e. 
^ ^ 0.5) and "relaxed" clusters in the turbulent energy - mass 
plane. As discussed in Section|4|the presence of "perturbed" clus- 
ters/groups introduces a bias in the properties of the overall sim- 
ulated cluster population. In this case the complete sample of our 
simulations would be more representative of rich environments and 
superclusters, with the smaller structures beeing more perturbed 
(and turbulent) than those in other environments. This would cause 
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a systematic flattening of tiie measured turbulent energy vs mass 
correlation. 



6 CONCLUSIONS 

Our sample of simulated galaxy clusters shows a well-constrained 
power-law scaling between the kinetic energy of ICM turbulent 
motions and the cluster mass. The slope of the scaling does not 
show any evident dependence on the resolution adopted to detect 
the pattern of turbulent motions, while the normalization shows a 
slight dependence on resolution. Rather independently of the clus- 
ter mass, the turbulence injected during the cluster life has an en- 
ergy budget of the order of 25% -35% of the thermal energy at 
redshift zero. If small perturbed clusters are included in the sam- 
ple this affects the scaling law of the turbulent energy making it 
slightly flatter than the thermal case; we expect that this should 
be the case in supercluster regions. Our scaling is in line with the 
semi-analytical findings of C&B05: although they used a com- 
pletely independent method, when plotted together the datapoint of 
the two approaches prove the scaling over more than two decades 
in cluster mass. Semi-analytical calculations give a simple physical 
explanation of the scaling laws in term of the PdV work done by 
the infalling subclusters through the main ones, and strenghten the 
physical nature of the turbulent field found in simulations. The in- 
clusion of cooling processes within our simulations is not expected 
to modify our conclusions, because the average cooling time for the 
large cluster regions considered here is longer than an Hubble time. 
Cooling may play an important role in innermost regions, where 
only a minor part of the turbulent energy is stored, however the in- 
clusion of cooling in simulations would also require the implemen- 
tation of feedback mechanisms - like galactic winds and bubble in- 
flation by AGNs - in order to prevent un-physical massive cooling 
flows. These processes should also introduce additional turbulence 
and further studies are required to understand how such processes 
might effect the reported correlations. 
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